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It is argued that the adjoint torelon loop, i.e. a Polyakov loop in the adjoint representation running in a spatial, 
rather than temporal, direction, is an observable which is sensitive to the presence of long color electric flux 
tubes at high temperatures. We show via lattice Monte Carlo simulations that this observable has a sharp peak 
at the deconfinement transition, remains much larger than the vacuum value for some range of T > T^, and falls 
below the vacuum value for T > 2Tc. This result suggests that long electric flux tubes may persist for a finite 
range of temperatures past the deconfinement transition, and at some stage disappear, presumably melting into 
a plasma of gluons. As a side remark, we point out that our results at T < Tc imply that the eigenvalues of 
ordinary Polyakov loop holonomies in the confinement phase have a slight tendency to attract rather than repel, 
which may be relevant to certain models of confinement. 



I. INTRODUCTION 

The thermal average of an operator (9 in a canonical en- 
semble is given by the standard formula of quantum statistical 
mechanics 



al — 



Y.„{n\0\n)e 



-E„/kT 



(1) 



where {|n)} is the set of energy eigenstates. In the case of a 
pure gauge theory, two questions naturally arise: First, at low 
temperatures, the thermal average is dominated by the ground 
state, so what is this ground state, exactly? The question could 
be answered, e.g. in the Schrodinger representation, by spec- 
ifying the Yang-Mills vacuum wavefunctional ^o[A] = (A|0). 
Secondly, at temperatures greater than the deconfinement tem- 
perature Tc, what are the relevant excited states that dominate 
the thermal average? If these relevant states |«) = Qn\0) at 
T > Tc can be constructed by acting on the ground state with 
ceitain gauge-invariant operators {2,,}, then what is the form 
of such operators? 

We are probably very close to having an answer to the first 
question, at least in D = 2 + 1 directions. It is likely, for the 
reasons given in refs. [1,2], that the Yang -Mills vacuum wave- 
functional has the form 
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or at least something quite close to that, where is the co- 
variant Laplacian in the adjoint representation, Aq is the lowest 
eigenvalue of —D^, and in is a mass parameter ' 



' An alternative possibility is that the kernel has instead the form 
1 l{\J —D^ — Ao + +'n), which reduces to the Karabali, Kim, Nair pro- 
posal [3] when evaluated on abelian configurations. The two kernels ai'e 
qualitatively very similar, and the corresponding wavefunctionals are diffi- 
cult to distinguish quantitatively in numerical tests, cf. [2]. 



The situation is much less clear at high temperatures, just 
past the deconfinement transition. All that can really be said is 
that in temporal gauge the relevant Q„ must be gauge-invariant 
operators, but this fact provides only modest guidance. Apart 
from such ultralocal operators as Tr[F^], gauge-invariant op- 
erators can be built out of combinations of the eigenvalues 
and eigenstates of gauge-covariant operators. For example, 
if is a Dirac operator, say, or a covariant Laplacian, and 
'Sl^n — X„(j)n is the corresponding eigenvalue equation, then 
any functional of the set of functions {(j)"{x)(j)l^(x)}, and the 
eigenvalues {A,,}, is gauge-invariant. The exponent in eq. (2) 
is another example. Gauge-invariant functionals can also be 
built from transformations to a physical gauge. Let g[x;A] be 
the transformation which takes a gauge-field A^. into Coulomb 
gauge. Then any functional G[g[x;A] oAi.(jc)] is also a gauge- 
invariant operator 

There are two types of operators, however, which may play 
a special role in pure Yang-Mills theory: the Wilson loops 
W{C), which create thin lines of electric flux, and the 't Hooft 
loops B{C), which create thin lines of magnetic flux known as 
center vortices. Unlike the previous examples, these operators 
do not require, as part of their definition, the solution of a dif- 
ferential equation over the full spatial volume. The low-lying 
energy eigenstates of strong-coupling lattice gauge theory are 
calculable, and can certainly be expressed in terms of Wilson 
loop operators acting on the vacuum state. There is also strong 
evidence that center vortices, created by the 't Hooft loop op- 
erator, play an important role in the confinement mechanism 
[4], including the generation of a spacelike string tension at 
high temperatures [5]. So one guess is that the relevant ex- 
cited states at high temperature are created by products and 
superpositions of these two special types of operators, i.e. 

|«) = ^{electric flux tube operators} 

X {magnetic flux tube operators} |0} . (3) 

Chernodub, Nakamura, and Zakharov [6] have presented evi- 
dence that a network of magnetic flux tubes in the deconfined 
phase is responsible for large contributions to relevant ther- 
modynamic quantities in the hot Yang-Mills equation of state. 

In this article we will be mainly interested in the electric 
flux tube content of the plasma state. At low temperatures. 
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well below Tc, the relevant Q„ will simply create low-lying 
glueball states, which are understood to be small closed elec- 
tric flux tubes. As T increases, longer electric flux tubes 
(highly excited glueball states), are less suppressed, and be- 
come increasingly relevant. It is a very old idea that flux tubes 
percolate at T = Tc', the flux tube free energy goes to zero 
while the energy per unit length remains finite. We are inter- 
ested in the question of what happens at T > T^. Do the color 
electric flux tubes simply vanish? Or do they persist for some 
range of high temperatures, before "melting" into a plasma of 
gluons? 



II. ADJOINT TORELONS AND ELECTRIC FLUX TUBES 

In this section we will explain why the thermal average of 
adjoint torelon lines, i.e. adjoint Polyakov loops in a spatial 
direction, might be relevant to the study of electric flux tubes 
at finite temperature. From the start we should confess that the 
argument we will present is entirely heuristic. Nevertheless, 
this argument serves to motivate our calculation of the adjoint 
torelon line, and it will suggest an interpretation of the results. 

We begin by considering Polyakov loops in the adjoint rep- 
resentation, of extension L, in the time direction. As is well 
known, adjoint Polyakov loops have a non-zero vacuum ex- 
pectation value even in the confined phase, because the color 
charge of a static adjoint source can be screened by gluons. 
The resulting state is known as a "gluelump," and it can be pic- 
tured in two ways. First, in some appropriate (e.g. Coulomb) 
gauge, a gluelump can simply be thought of as a bound state 
consisting of the static source and a single constituent gluon. 
Fig. 1(a). Secondly, just as a glueball can be pictured as a 
closed, fluctuating string of electric flux, a gluelump can be 
thought of as a closed electric string which begins and ends at 
the static source. Fig. 1(b). 

The evidence that a closed electric flux tube can be de- 
scribed by a fluctuating string has become quite compelling 
in recent years, for both 2h-1 and 3+1 dimensions. In particu- 
lar, lattice Monte Carlo simulations have found that the low- 
lying spectrum of closed flux tube states, which are closed by 
lattice periodicity, is described to a surprising degree of ac- 
curacy by the simple Nambu-Goto action. The ground state 
energy in particular agrees with the Casimir energy of the 
Nambu string (associated with the Luscher term) down to tube 
lengths where one might think that a simple string picture of 
the flux tube should not apply. Up-to-date results can be found 
in refs. [7, 8]. We will therefore be guided by the picture of 
a glueball as a closed, fluctuating string, and assume that this 
picture also applies to gluelumps, which would then resemble 
strings attached to static DO-branes. A gluelump has some fi- 
nite extension which, like that of most low-lying hadrons, is 
probably on the order of one fermi. This rough estimate is 
supported, to some extent, by the fact that the string between 
two static adjoint sources breaks into two gluelumps when the 
sources reach a separation of 1.25 fm, as determined from lat- 
tice Monte Carlo studies [9]. The length of the "broken" string 
attached to one of the sources, at the breaking point, serves as 
a crude upper bound for the radius of the gluelump. 



Now consider what will happen as the extension L- of the 
periodic volume in the spatial z-direction is reduced to or be- 
low the extension of the gluelump. In the constituent gluon 
picture, the bound gluon can take advantage of the periodicity 
in the z-direction by having zero momentum in this direction; 
this reduces the kinetic energy of the bound gluon without 
increasing color charge separation beyond the normal exten- 
sion of the gluelump. The ground state gluelump wavefunc- 
tion is then completely delocalized in the z-direction. Like- 
wise, in the string picture, the string can wind around the pe- 
riodic z-direction, as indicated in Fig. 1(c), without incurring 
a large cost in energy. In the string picture, as becomes 
very small, the string partition function is dominated by states 
with multiple windings, and the contribution from states with 
zero winding number is negligible by comparison. In the case 
of the gluelump, where the ends of the string are pinned to a 
source and the center of mass cannot propagate freely in the 
z-direction, the winding mode represents an extra degree of 
freedom, strongly suppressed at large L,, which opens up at 
small L,. 

In string theory, of course, the L- limit of a DO-brane 
is equivalent, by T-duality, to the L- — > oo limit for a Dl-brane. 
In the latter case, summation over the center-of-mass momen- 
tum in the z-direction contributes a factor of to the ther- 
mal partition function (see, e.g., Zwiebach [10]), so the Dl- 
brane partition function diverges linearly with L-. Likewise, 
the sum over winding modes of the DO-brane leads to a di- 
vergence as 0. Of course, the expectation value of an 
adjoint Polyakov line at finite temperature is not necessarily 
the same thing as the thermal partition function of string end- 
ing on a DO-brane. Nevertheless, since the low-temperature 
flux tube spectrum is so well described by the Nambu action, 
we are led to expect that string winding modes will enhance 
the expectation value of the adjoint Polyakov loop, as L~ de- 
creases. 

Next let us interchange the z and t labels, so that the "short" 
direction is now the time direction, and the Polyakov loop runs 
along the large spatial z-direction. Polyakov loops which run 
in a spatial direction are sometimes referred to as "torelon s," 
and we will adopt that terminology here. The expectation 
value of an adjoint Polyakov loop in a 4-volume with one short 
space direction can be reinterpreted as the thermal expecta- 
tion value of an adjoint torelon loop at finite temperature. The 
string pictures of the gluelump state. Figs. 1(b) and 1(c), now 
take on a somewhat different meaning. 

Before going into that, it is useful to first note that physical 
states in an SU(A^) pure gauge theory can be characterized by 
their transformation properties with respect to an underlying 
global center symmetry defined at any fixed time. In a lattice 
formulation, in a f = timeslice, consider the global transfor- 
mation 

U-Sx,y,z = 0) ^ e'2'^/"f/,(x,>',z = 0) , (4) 

with all other links at z 7^ unchanged. A physical state '^\U] 
has A^-ality k if, under the global transformation (4), 

^>[U] e'2'^*/'v^[[/] . (5) 

An operator which changes the A^-ality of a state by M = 1 is 



static source static source 
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FIG. 1. Descriptions of the gluelump state: (a) The gluelump as a bound state of a static source and a constituent gluon in Coulomb gauge. 
Dashed lines indicate the color Coulomb interaction. The red/blue shadings correspond to the adjoint representation as a two-index object, (b) 
The gluelump as a fluctuating electric string pinned to the static source, (c) For spatial extension L- comparable or smaller than the gluelump 
size, the electric string will tend to wind around the periodic volume. 
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FIG. 2. Schematic diagram of string propagation in a space-slice 
at fixed z- (a) Creation at f = of an electric string-antistring pair, 
which subsequently annihilate, (b) An electric string created at f = 
which propagates through the periodic lattice, and is destroyed at 
/ = 0. 



t=Lt/2 





t = 



(a) 



(b) 



FIG. 3. Lattice strong-coupling diagrams corresponding to (a) Fig. 
2(a); and (b) Fig. 2(b). 



the torelon operator in the fundamental representation 



^ = TrPexp 
or in lattice regularization 

^ = Tr 
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while the conjugate operator changes the A^-ality by Ak - 
— 1 . Any superposition of Wilson lines 



^fl(C)TrPexp 
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which are closed by lattice periodicity in the z-direction, with 
winding number equal to one, also change the A^-ality of states 
by one unit. A thick flux tube which winds once through the 
periodic lattice is created by operators of this type, acting on 
the vacuum. The vacuum state itself, assuming ordinary pe- 
riodic boundary conditions, has A^-ality zero, and in fact the 
lattice version of the vacuum wavefunctional (2) is explicitly 
invariant under the transformation (4). 

Because ^ and take states from one A^-ality sector to 
another, we can regard them as creating/destroying long elec- 
tric strings which wind through the periodic lattice in the z- 
direction. Operator 3^ creates a string or destroys an "anti- 
string," i.e. a string with the opposite orientation of color elec- 
tric flux, while creates an antistring or destroys a string. 
The adjoint torelon loop is simply a Polyakov loop which runs 
in the z-direction, in the adjoint representation of the gauge 
group, i.e. 



(9) 



Now the string diagrams for gluelumps. Figs. 1(b) and 1(c), 
under the z ^ f interchange, can be reinterpreted as follows; 
At low temperatures, with L, much greater than the gluelump 
extension. Fig. 1(b) can be viewed as the creation of a string- 
antistring pair, which propagate for a finite time and subse- 
quently annihilate. This process is indicated schematically in 
Fig. 2(a), which shows a slice of the volume at fixed z; i.e. 
the adjoint torelon line runs into the page. The correspond- 
ing diagram in lattice strong-coupling perturbation theory is 
shown in Fig. 3(a); it involves a sheet of plaquettes wrapping 
around the adjoint line to form a tube. It is contributions of 
this type that contribute to the vacuum expectation value of 
the adjoint torelon loop. At high temperatures, with less 
than the typical gluelump extension of roughly one fermi, the 
winding string in Fig. 1(c) becomes relevant. This figure can 
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FIG. 4. Thermal average of the adjoint torelon hne (^A)thential on a IS-' x L, lattice, vs. temperature ratio T /Tc. The deconfinement temperature 
at j3 = 2. 5115 corresponds to L, = 8 lattice spacings in the time direction, and therefore T /Tc = S/Lt. The dashed line is inserted to guide the 
eye. 




be reinterpreted, as in Fig. 2(b), as a string worldsheet which 
wraps around the periodic time direction, passing through the 
z-axis at f = 0. The corresponding lattice strong-coupling dia- 
gram is shown in Fig. 3(b); in this case the sheet of plaquettes 
wraps around the periodic time direction. As explained in, 
e.g., ref. [11], a particle worldline which wraps around the pe- 
riodic volume in the time direction can be interpreted as a real 
particle in Minkowski space, rather than a virtual particle in 
Euclidean space, which contributes to the thermal ensemble. 
Likewise, the wrapped worldsheet represents contributions to 
{^A)rhermai due to States containing at least one long electric 
string, closed by periodicity in the z-direction. 

The identification of hyperplanes at f = ±Lr/2 corresponds, 
in the canonical ensemble, to the sum over states in (1). We 
have already argued that the degree of freedom associated 
with the string winding mode, which opens up at small com- 
pactified dimension, ought to enhance the expectation value 
of an adjoint loop. But since a winding number ^ > corre- 
sponds, in the thermal field theory, to states of A^-ality k mod 
A^, such processes can only contribute if states containing a 
minimum of k long flux tubes are present in the thermal en- 
semble at high temperatures. Thus, the greater the electric flux 
tube content in the thermal ensemble, the greater the contribu- 
tion from the wrapped worldsheets, and the larger the value 
of {^A)thermai- Conversely, if long electric flux tube states are 
for some reason almost absent in the thermal ensemble, then 
only worldsheets with zero winding number and associated A^- 
ality A: = can contribute, as in the T = case, and all other 
potential contributions are suppressed. 

The qualitative argument just presented is our reason 
for thinking that any significant enhancement of an adjoint 
torelon line, beyond its T = value, may indicate the pres- 
ence of states containing long electric flux tubes in the ther- 
mal ensemble. With this motivation, we have computed the 
expectation value of the adjoint torelon line in SU(2) lattice 



gauge theory, using the standard Wilson action at )3 = 2.51 15 
on lattices of volume 15-' x L,. The extension L, in the time- 
like direction runs from 15 down to 2 lattice spacings, and 
the torelon line, running in the z-direction, is always 15 lat- 
tice spacings long. The critical temperature at j3 = 2.5115 
corresponds to 8 lattice spacings in the time direction [12], so 
we are sampling temperatures from a low of j^Tc to a high at 
ATc- The lattice spacing at )3 = 2.5115 is a = 0.08 fm, so if 
the extension of the gluelump is about one fermi, we would 
expect that winding contributions would become significant 
at around 12 lattice spacings {T = ^Tc). In this calculation 
we have made use of the Liischer-Weisz noise reduction tech- 
nique [13] with one level of sublattices, each of three lattice 
spacings extension in the z-direction. 

Our results are displayed in Fig. 4. Although the mag- 
nitudes of the adjoint torelon expectation values are always 
very small compared to (9(1), there is a dramatic structure 
seen in the relative values, which show a sharp peak precisely 
at the deconfinement transition. According to our previous 
reasoning, this peak could be attributed to the presence of 
long electric strings for temperatures at or near the decon- 
finement transition. Two points are worth noting. First, the 
adjoint loop does not immediately drop to (or below) the vac- 
uum value at r > Tc- For temperatures a little beyond Tc, 
the thermal average is still much greater than the vacuum ex- 
pectation value. On the other hand, as T increases beyond 
TTc, {^a) thermal f^Hs to valucs which are very much smaller 
than the low-temperature value. Our interpretation of these 
results is that states with long electric flux tubes become in- 
creasingly important as the temperature increases, up to the 
deconfinement transition temperature. Beyond that transition, 
long electric flux tubes persist for some range of temperatures, 
until, at high enough temperatures, the relevant excited states 
no longer contain such configurations, and the flux tubes pre- 
sumably "melt" into the surrounding plasma. 
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III. DO POLYAKOV LOOP EIGENVALUES ATTRACT, OR 
REPEL? 

Our calculation of the adjoint torelon line at the lowest 
temperature T ~ j^Tc, obtained on a hypercubic 15"* lattice 
volume, is also (again interchanging the z and f-directions) 
a reasonable estimate of the expectation value of an adjoint 
Polyakov loop at this (comparatively) low temperature.^ This 
data point allows us to address quite a different question, 
namely: In a pure Yang-Mills theory, in the confinement 
phase, do the eigenvalues of a fundamental representation 
Polyakov loop tend to attract, or tend to repel? 

Obviously the endpoint of attraction is where the eigenval- 
ues are identical, and this must be a center element, while 
the endpoint of repulsion is a situation where the trace of the 
Polyakov loop is zero. An example of the former is Zn lat- 
tice gauge theory, where Polyakov loops are precisely center 
elements and the VEV vanishes, in the confined phase, due 
to quantum fluctuations among these center elements. This is 
also the case for center-projected SU(A^) lattices in maximal 
center gauge. The situation is similar to the disordered phase 
of a Zn spin system where there are islands of the same cen- 
ter element, but upon averaging over the full volume, the ele- 
ments cancel. In the ordered phase, the lattice is overwhelm- 
ingly one center element or another The opposite extreme is 
the dilute dyon gas advocated by Diakonov and Petrov [15]. 
In that picture, in the confined phase, the trace of almost any 
Polyakov loop (away from the middle of the dyon) is zero in 
each configuration, which is the limiting case of eigenvalue 
repulsion. Deconfinement occurs when there is a switch (at 
the minimum of the free energy) from eigenvalue repulsion to 
eigenvalue attraction. 

Of course, these pictures are extremes. In reality the dis- 
tribution of Polyakov eigenvalues is close to random in the 
confined phase, i.e. close to the Haar distribution.^ The same 
can be said for the eigenvalue distribution of large Wilson 
loop holonomies. But there is, nevertheless, a lot of informa- 
tion contained in the eigenvalue distribution of Wilson loop 
holonomies. The tiny deviation from the Haar measure in that 
case encodes information about the string tension, the Liischer 
term, Casimir scaUng, and color screening. For this reason, 
we may expect the small deviation from the Haar measure for 
Polyakov loops to be significant as well. If the underlying 
confining configurations tend to have Polyakov eigenvalues 
coincide, as in a center vortex mechanism, or repel, as in a 
dilute dyon gas, then these tendencies might be expected to 
show up in the small deviation of the eigenvalue distribution 
away from the Haar measure. 



Let 



f/(L) =Pexp 



A previous calculation by Gupta et al. [14] of the adjoint Polyakov line 
in SU(3) pure gauge theory was carried out in the confinement phase for 
temperatures quite close to the transition temperature, T/Tc > 0.907. A 
positive value was obtained, but for our purposes it is desirable to check 
this result at a much lower temperature, well below the transition point, 
where high temperature effects are hopefully negligible. 
' In fact, even if center symmetry were explicitly broken by matter fields in 
the fundamental representation, the eigenvalue distribution would still be 
nearly random at low temperatures, just due to non-confining fluctuations. 



dtAo{0,t) 



(10) 



be a Polyakov loop holonomy, at spatial coordinates x taken 
to be X = 0, and let P = Tr[f/(L)]. If the eigenvalues of U{L) 
attract, then the probability distribution for U (L) on the group 
manifold should be slightly peaked around center elements; 
i.e. P = zN, where z ^ Zn for the SU(A^) gauge group. Con- 
versely, if the eigenvalues tend to repel, then the distribution 
will be peaked around configurations for which P = 0. 

The probability density p{g) for a Polyakov holonomy is 
given by 



p{g) = {5[g-U{L)]). 



(11) 



where the Haar measure distribution corresponds to p (^) = L 
Now make a character expansion of the delta function, and at 
this point specialize to the SU(2) gauge group, so that 

5[g-UiL)]^ £ Xj[U{L)]Xj[g]. (12) 
j=o,l,i,... 

Then, in the confined phase 

{d[g~um^ E {xjwmxjM- (13) 

j=0,l,2,3... 

Note that the half-integer representations are dropped in the 
confined phase, on grounds that {Xj[U{L)]) vanishes, for j = 
half-integer, if center symmetry is unbroken. But, apart from 
the trivial j = representation, Polyakov loop VEVs are very 
small even for integer j, and are very rapidly suppressed as 
j increases. It is therefore an excellent approximation to the 
sum to keep only the leading contributions, i.e. 

p{g) = {8[g^U{L)]) ^l + {TiA[UiL)])TiA[g] , (14) 

where Tr^ [g] = Xi [g] is just the trace in the adjoint represen- 
tation. 

The crucial question is now whether {Tia[U{L)]) is posi- 
tive, or negative. If it is positive, then the Polyakov loop prob- 
ability density is largest when g is a center element (= ±1 for 
SU(2)), and we can conclude that the eigenvalues tend to at- 
tract. Conversely, if the trace of the adjoint Polyakov loop is 
negative, then the probability density is largest at Tr^ [g] = — 
which is obtained when the trace in the fundamental represen- 
tation is zero. In that case, the eigenvalues tend to repel. 

We have seen, from the first data point (obtained on a 15^ 
lattice) shown in Fig. 4, that the trace of an adjoint Polyakov 
loop at low temperatures is positive. We conclude from this 
data that the eigenvalues of a Polyakov loop tend to very 
slightly attract, rather than repel, in the confinement phase. 



IV. CONCLUSIONS 

If the qualitative argument in section II is correct, then the 
peak we have seen in the adjoint torelon loop at the deconfine- 
ment transition is an indication that long electric flux tubes 
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survive, in the thermal ensemble, somewhat beyond the de- 
confinement transition, up to T « 2T^.. For a true quark-gluon 
plasma created via heavy ion collisions, this would suggest 
that some hadronic bound states, such as closed flux tubes 
(i.e. glueball states) and certain quark bound states, may sur- 
vive a little beyond Tc- This conclusion ties in very well with 
the discussion of Shuryak in ref. [16], who argues that certain 
aspects of recent RHIC data could be explained by the exis- 
tence of metastable color electric flux tubes in the quark gluon 
plasma, at temperatures somewhat beyond T^. Shuryak and 
co-workers have also argued for the existence of flux tubes 
past the transition based on a dual superconductor model of 
confinement, cf. [17, 18]. As we will show in a subsequent 
article [19], the survival of electric flux tubes at T > Tc, in 
the context of the gluon chain model [20], may also be con- 
nected with the persistence (indeed, the increase) of the color 
Coulomb string tension in the deconfined phase. 



We have also seen that the mere fact of positivity, of ad- 
joint torelon loops at low temperatures, implies that the eigen- 
values of fundamental representation Polyakov loops tend to 
slightly attract, rather than repel, in the confinement phase. 
If this small deviation from the Haar measure reflects the na- 
ture of the underlying confinement mechanism, then confin- 
ing configurations associated with center elements, i.e. center 
vortices, would seem to be favored. 
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